Molecular Correlations in a Supercooled Liquid 
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We present static and dynamic properties of molecular correlation functions Si mn> ii m i n i (q,t) in a 
simulated supercooled liquid of water molecules, as a preliminary effort in the direction of solving 
the molecular mode coupling theory (MMCT) equations for supercooled molecular liquids. The 
temperature and time dependence of various molecular correlation functions, calculated from 250 
ns long molecular dynamics simulations, show the characteristic patterns predicted by MMCT and 
shed light on the driving mechanism responsible for the slowing down of the molecular dynamics. We 
also discuss the symmetry properties of the molecular correlation functions which can be predicted 
on the basis of the C^-symmetry of the molecule. The analysis of the MD-results for the static 
correlators Si mn jt m i n t(q) reveals that additional relationships between correlators with different 
signs of n and n' exist. We prove that for molecules with C rl ,-symmetry this unexpected result 
becomes exact at least for high temperatures. 

PACS numbers: 64.70.Pf, 6f.25.Em, 6f.20.Ja 



I. INTRODUCTION 

In recent years, significant progress has been made 
in the understanding of the slow dynamics in super- 
cooled liquids @. Theoretical [|-|, experimental (see 
e.g. |P-|To|| ) and simulation efforts (see e.g. 11-14Q have 
highlighted the role played by the ideal-glass transition 
temperature T c , first predicted by mode-coupling the- 
ory (MCT) for simple liquids fll5| , [l6| and in a schematic 
model [ p"7| , and identified the different dynamical mech- 
anisms above and below this temperature. This work 
has clarified the strong interplay between liquid struc- 
ture and liquid dynamics above T c as well as the uni- 
versal aspect of the decorrelation process — which close 
to T c is predicted to become independent of the corre- 
lation function as well as the q vector that is probed 
fLsf . For colloidal systems MCT has been tested both 
on a semi-quantitative level ||, i.e. the validity of the 
scaling laws was investigated, and also on a quantita- 
tive level where it has been demonstrated that the time- 
dependent density correlator for a liquid of hard spheres 
(which is a good model for neutral colloids) obtained 
from MCT describes the corresponding experimental re- 
sult over 3 decades in time by using only one fit param- 
eter (for the timescale) |T^]2C[] . All the other experimen- 
tal tests were semi-quantitative and were exclusively re- 
stricted to molecular systems, like OTP, Salol, Glycerol, 
etc. (for more details see refs. [|]||], ]7|~ |l0| , ^l| and ref- 
erences therein). On a first glance the reasonably good 
agreement between the predictions of MCT for simple liq- 
uids with many of the experimental data for the molec- 
ular glass forming liquids seems to be surprising, since 
the orientational degrees of freedom do not appear in 
the original version of MCT [R~8| . To do a quantitative 



comparison between theory and experiment or simula- 
tion and especially to describe the molecular correlations 
(including orientational degrees of freedom) it is neces- 
sary to extend MCT to molecular liquids. In particular, 
a molecular mode coupling theory (MMCT) allows to 
study the role of the coupling between the translational 
degrees of freedom and orientational degrees of freedom. 
Such an extension has recently been performed for a sin- 
gle dumbbell molecule in an isotropic liquid [p2[ and for a 
liquid of diatomic molecules |p3|-p5|. Reference and 
p5| even treat general molecules. 

MMCT is also conceptually based on the hypothesis 
that the liquid structure is controlling the long-time dy- 
namical evolution of the system. Structural information 
is used as input in the theory via generalization of the 
density-density structure factor including the angular de- 
grees of freedom. Such quantities, although they are 
difficult to determine experimentally, can be evaluated 
from molecular dynamics trajectories and used to test 
the quality of the recently proposed MMCT. For the case 
of a liquid of diatomic molecules, the molecular structure 
factors have been evaluated |2(| and a first quantitative 
test of MMCT has been presented in Ref. (2^j2|. Fur- 
thermore, the temperature dependence of the molecular 
correlation functions can give hints regarding the leading 
mechanism for the slowing down of the dynamics. 

The asymptotic predictions of MCT continue to be 
valid within the MMCT scheme. Also in MMCT |22]j23| 
i) there exists a /3-relaxation regime where the factor- 
ization of the time- or frequency-dependence of the cor- 
relators or susceptibilities from the space- and angular 
dependence holds generically and the time-dependence 
is given by the /3-correlator which fulfills the first scaling 
law fl8f| and ii) there exists an a-relaxation regime in 
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which the second scaling law Q holds, i.e. the 
time-temperature superposition principle is fulfilled [^9) . 
These two results underline the universality of the ideal 
structural glass transition based on a bifurcation scenario 
described by a fold-singularity |Q. However, the cal- 
culation of the various exponents (entering the scaling 
laws), the calculation of the critical nonergodicity param- 
eter, the critical amplitudes, the transition temperature 
T c and particularly of the molecular correlation functions 
themselves require the solution of the MMCT equations. 

In this Article we present the temperature and q- vector 
dependence of the generalized structure factors required 
by the MMCT for the case of SPC/E potential @. This 
potential, which describes the molecule as a rigid planar 
body and models the pair interactions as a sum of elec- 
trostatic and Lennard Jones terms, has been studied in 
detail and has been shown to reproduce qualitatively the 
characteristic properties of liquid water — a liquid where 
the slowing down of the dynamics on cooling is not re- 
lated to packing constraints but to the formation of a 
tetrahedral network of highly directional hydrogen bonds 
pl|p2| . Previous numerical studies on the SPC/E sys- 
tem, focused on the center of mass self |}3| and collective 
properties fl34j , have shown on a semi-quantitative level 
that the center of mass dynamics is well described by 
MCT, with an estimated critical temperature of about 
T c = 200 ± 3K. It has been shown that orientational 
degrees of freedom are crucial in strongly enhancing the 
slowing down of the dynamical processes on supercool- 
ing. This has led to the introduction of a semi-schematic 
MCT model ]|^] where the coupling to the orientational 
degrees of freedom is accounted for phenomenologically 
by the introduction of a parameter \R > 1 by which the 
coupling of the center of mass-density modes is enhanced. 
This semischematic approach, recently tested also for the 
case of a supercooled liquid of diatomic molecules [j36| , 
focuses only on the center of mass translational degrees 
of freedom and is by construction unable to describe the 
time evolution of the orientational degrees of freedom and 
to indicate which orientational degrees of freedom are re- 
sponsible of the slowing down of the center of mass and 
orientational dynamics. A MMCT description is required 
to fully describe these important dynamical aspects. 

The generalized structure factors presented in this Ar- 
ticle are a first step in the direction of solving the MMCT 
for a molecular system. In section [n| we motivate the in- 
troduction of an infinite-dimensional correlation matrix 
Si mn ,i'm'n'(q,t) where I = 0,1,2, ... ; -I < m < I , —l< 
n < I and corresponding relations for the primed quan- 
tities. As already mentioned above, MMCT requires the 
static correlators Si mni i> m ' n '(q) as an input. In Sec. Ill 
we discuss the relations between the Si mni i> m > n >(q,t) cor- 
relation functions arising from the symmetries character- 
istic of isotropic liquids and the symmetries characteris- 
tic of C rv molecules. These relations reduce drastically 
the number of independent Si mni i> m ' n '(q,t) correlation 



functions and support the feasibility of a full MMCT 
calculation for this class of liquids. Finally, in Sec. IV 
we present the time evolution of the generalized angular 
correlation functions for the SPC/E case and interpret 
their behavior in the general MCT framework. The qual- 
itative and semi-quantitative agreement between numer- 
ical data and asymptotic MCT theoretical predictions 
strongly suggests to perform a full MMCT comparison. 



II. MOLECULAR CORRELATION FUNCTIONS 

To describe a molecular liquid it is necessary to in- 
troduce, besides the information on the position of the 
molecule's center of mass, the information on the orien- 
tation of the molecules. The microscopic density, defined 
for simple (atomic) liquids (3^] as p(x,t) ~ 5Z 7 =i — 



j(t)), is generalized to 



N 

£ 



p(x, n, t) = V 6(x - xj {t))s(n - rij (*)) (i) 



where the sum runs over the N molecules of the liquid, 
Xj(t) and flj(t) = (4>j(t),9j(t), Xj(t)) are respectively the 
position of the center of mass and the Euler angles of the 
jth molecule at time t. 

Any function f(x,Q) can be expanded with respect 
to plane-waves and to generalized spherical harmonics 
D l mn (n) = e-*«* d l mn {9) e-*"x ( see Ref. §§) as 



/(^)^EEH/(2i+i)- 

q Irnn 

■ flmn(q) e~^ S D l mn (Q) 



(2) 



where the coefficients fi m n{q) are given by 



flmn(q) = i l (2l + 1)> / d 3 x / dSlf(x,tl) e iq * D^ n (n) 



(3) 

Here it is I > 0, —l<m<l, —l<n< I. Application of 
@ to p(x, J7, t) yields the tensorial one-particle density: 

N 

p lmn (q,t) = i\2l + l)' £ e W) D^Mt)). (4) 

The prefactor in eq.(^) is chosen for technical conve- 
nience, e.g. the factor i l makes the molecular correla- 
tors real for n — n' . As will be shown below, this choice 
produces real correlators also for n and n' ^ at high 
temperature. Now we can introduce the time-dependent 
molecular correlation functions 

Slmn,l>m'n> = {P* m n (& *) Pl'm'n' (q, 0)) (5) 
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where the angular brackets denote the canonical average 
over the initial point in phase space. Substitution of (^) 
into (H) yields 

Slmn,l'm'n'(q,t) = S\2n,l'm>n> *) + ^Jn.l'roW®*) (*>) 

with the distinct part 



and the se// part 



j 

■ (e-*iW-^)C(%(i)) ^n'(%)>- (8) 

The reader should note that these correlation functions 
involve both translational degrees of freedom and ori- 
entational degrees of freedom. They form a complete 
set for any two point correlation function of an arbitrary 
molecular liquid and they are also the main entities en- 
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tering MMCT 



Specialization to n 



yields 



the corresponding correlators used in Ref. |23| for lin- 
ear molecules. Although their introduction is enforced 
by theoretical reasons since homogeneity and isotropy of 
the liquid is accounted for directly by the Fourier trans- 
formation to the g-space and their tensorial nature, only 
very few of them can be measured experimentally. For 
instance St mn> i mn (q = 0,i) and S^ n2mn {q = 0,t) can 
be obtained from dielectric measurement and NMR, re- 
spectively. Information on the center of mass-correlator 
Sboo,ooo(<f, t) follows from light scattering, provided that 
the contribution of the orientational correlators can be 
neglected [[39|. The neutron scattering cross section is a 
linear superposition of all of these correlators (see e.g. 
p0|) from which information on the individual correla- 
tors could be obtained by choosing different scattering 
lengths of the atomic units. On the other hand, it is 
a great advantage of a MD-simulation, which really de- 
termines the trajectories {xj(t), fij(i)}, that these cor- 
relation functions can be calculated. This of course can 
only be done for I and I' smaller than a cut-off value l co 
which will in the analysis of the SPC/E data be chosen 

tlrS l QQ 2 . 



III. SYMMETRY PROPERTIES 

In this section we will discuss the general properties of 
the correlators (^),(0) and (ph which follow from symme- 
try. Similar discussions have been done for an expansion 
into rotational invariants |Q and in real space . The 
properties presented in the following will be of great im- 
portance for the discussion of the results of the simulation 



and especially to reduce the effort for the solution of the 
MMCT equations which we will present in a subsequent 
paper. Concerning symmetry, we have to distinguish be- 
tween global and local symmetries where the latter are 
related to the geometry of a single (rigid) molecule. 



A. Global Symmetry 

The global symmetries arise from the invariance of the 
molecular interactions under the simultaneous transla- 
tion or rotation of all molecules, provided that the ex- 
ternal potentials arc zero. Similarly the absence of a 
time-dependent external force implies time translational 
and time reversal symmetry which are the same as for 
simple liquids and therefore will not be discussed here. 
The invariance under translations in space has been al- 
ready accounted for by the transformation into (f-space, 
i.e. what remains is the discussion of the transformation 
of the molecular correlators under rotations and reflec- 
tions Q. 

Let us start with the proper transformations 
R(a,(3, 7) G SO(3). Immediately from the transforma- 
tion of D l mn (Q) under R follows the transformation 
law for the molecular correlators: 

Si mn ,l'm'n'(q',t) = E D m"m( R ) D m"'m'{ R ) ' 
m" ,m"' 

■ Sl m "n t li m i" n i(q,t) (9) 

where q' = Rq. The choice of the g-frame, i.e. the labo- 
ratory frame of reference oriented such that q lies along 
the z-axis, q = qo = (0,0, g), leads to a simplification 
of the molecular correlators. For a rotation R z (a) = 
R(a,0, 0) around the z-axis for which q' = R z q}> = qo, 
eq.(||) implies: 



S, 



mn.V m'n 



(90,*) 



-i(m-m')a S lmn>Vm , n ,fo,t) (10) 



where Di m ,( R z (a)) = e 



38 1 has been used. 



Since ( |10[ ) is valid for all a, the correlators in g-frame 
must be diagonal in m and m': 



(qo,t) = S 

In.l'n' 

fat) 8 



(ii) 



where q= \q \. Since (|9|) is also true for S\^ n>l , m , n ,{q,t) 

and s\^ inVm , n ,{q,t) the m-diagonality also holds for 
these quantities. A further relationship follows from (||) 
by choosing a rotation Ry(ir) = R{0, tt, 0) by 7r around 
the y-axis: 

£>lmn,l'm'n' ( <70j *) = ( — 1) Slmn.l'm'n' (ffO) *)■ 

(12) 

Here we considered that Ry(%)qo — —qo, D l mn (Ry{%)) = 
(— \) l+m S mn [B8| and m (or n) denotes — m (or — n). 
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Next we investigate the inversion P for which 
D l mn (Pn) = (-1)'+" D\ n JSl) @. Then © yields 

Slmn,l'm'n'( — q,t) = (— 1)'+"+' + n Sl m n,l' m' n> (<7> *) • (13) 

which is valid for arbitrary (f. 

A further useful relation is obtained taking the com- 
plex conjugate of S7 m „ i' m 'n' (<f, i) and applying the equal- 
ity D£ n (fi) = (-l)™+»Dj„ n (fi) @: 



for i/ an integer and /irced j implies that 



2tt 



(21) 



[*S/ mn ,/'m'n'(^)]* = (~ 1) 



+m+m' -\-n-\-n 



Sim. n,l'm'n' ( (?j ^) 

(14) 



Combining eqs.(14) respectively with eq.Jl^) and eq.jl^) 
we have the final relations: 

[Slmn,l'm'n'(<l,t)]* = (~ l) m+m SlmnMm'n'(q,t) (15) 

[SWm'n'&t)]* = (-l) n+n 'S im „^ m , 2 ,(£i) (16) 
or in the q-frame, due to the diagonality in m: 

ISJZMfaW = St Vn ,(9o,t) (17) 



[s%,vn'(qo,t)]* = (-i) n+n siivAM (18) 

We will show in the following section that, in the case of 
molecules with C% v symmetry, like water, n + n' must be 
even and, thus, it holds the symmetry 

[SZM^t)]* = S^ lln ,(q ,t) (19) 



B. Local Symmetry 

In case that a molecule has a point symmetry we can 
derive additional identities for Si mn> i' m > n >(q,t). Since we 
have water in mind which has a C2„-symmetry we will 
discuss molecules with C rv symmetry (see e.g. To 
avoid confusion with the n-index of the rotation matri- 
ces D l mn we deviate from the conventional notation C nv . 
This symmetry means that the molecule possesses a r— 
fold rotational axis and r planes of reflection symmetry 
which contain the rotational axis. It is obvious that ap- 
plication of the local symmetry operations to any single 
molecule must leave the interaction invariant. Let us 
begin with the r-fold rotational symmetry. Without re- 
striction of generality we choose the body-fixed z-axes 
along the r-fold symmetry axis. In that case the r-fold 
symmetry only affects the third Euler angle x- The trans- 
formation 



Xj 



2tt 

Xj + v — 

r 



(20) 



Since D l mn (Qj(t)) = e~ im *^ d l mn (Mt)) e -™>0«) || 
it follows with the separate use of (20) and (21) to the 
distinct part (|^) : 

e iu ^ n = 1 and e"^™' = 1 (22) 

for all integers v. This restricts n and n' to integer mul- 
tiples of r, i.e. 

q {d) ,..{^0 forn,n'G{0,±r,±2r,...} 

°lmn,l'm'n'\.Q> > | Q otherwise 

(23) 

For the self part (eq.(||)) we have to use ( |20| ) and ( |2l| ) 
simultaneously, leading to: 



^-(n-n') _ ^ 



(24) 



for all integers v, which is fulfilled for n' = n(mod r), i.e. 

w r „- n = /^° for(n-n') 6 {0, ±r, ±2r, ...} 
J imn,i'm'n'\q^ ) | q otherwise 

(25) 

According to these results the total correlator 
Simn,i'm'n'(q, t) (eq.^)) reduces to its self part for n or n' 
not equal to an integer multiple of r. For the case of wa- 
ter (r=2), eqs.(|24|) requires n + n' to be even, a condition 
which simplifies cq.(|l8|) further. 

What remains to be discussed is the role of the reflec- 
tion symmetry. Here we restrict ourselves to the static 
correlators. It will be shown in Appendix |X] that for their 
distinct part it implies: 



(d) 



T — > oo 



(26) 



i.e. it does not depend on the sign of n and n'. 

We close this section by a discussion of the implica- 
tions following from the results in both subsections. It 
will be crucial that we will consider the correlators in the 
q-frame. The correlators in an arbitrary reference frame 
are easily obtained from eq. ([)]). Since the identities de- 



rived in subsection [ill A| also hold for S^ n Vm , n , (q, t) and 



S, 



(d) 



■in, I'm' n 



(q,t) we get from eqs.(0)-mj) with (E3J) and 



(Eq) for static correlators in case of high temperatures: 



{d)m 



-y(d)m 



(27) 



i.e. the static distinct part is real and does not depend 
on the sign of m. Since the static self part fulfills: 



S. 



(s)m 
In.l'n 



(q) = S U > S„ 



(28) 
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it follows from (g), (g6|) and ( p8| ) for high temperatures 
that 

,(d)r 



where 



S I|nl,i'|n'|(8) ' 71 f 9Q ^ 



^ini'n'il) ^ S approximately determined by 
,| (g) and in additioi 
on the sign of ra, due to (E 



^ilnlT'ln'l (?) an< ^ * n Edition it is real and does not depend 



IV. RESULTS 

In this section we present a detailed analysis of the 
static and dynamic rotational correlators as calculated 
from a Molecular Dynamics (MD) simulation of a molec- 
ular network- forming liquid. The system is a liquid of 216 
rigid molecules whose geometric parameters are chosen in 
such a way to mimic water molecules. The intermolecu- 
lar interactions are described by the SPC/E (Single Point 
Charge Extended) potential jIo| which has been shown 
to be able to describe most of the thermodynamic prop- 
erties of supercooled liquid water. We will not report 
here the details of the simulation which can be found in 
refs. US 

The low T simulations have been run for more than 
50 s integration timesteps, corresponding to 50 ns. Each 
of the simulation runs requests several months of CPU 
time on one 400 MHz alpha-processor for equilibration 
and more than one year of computer time to generate the 
configuration ensemble studied. In this situation, equi- 
librium runs for much larger systems at low temperature 
are very time consuming. Since periodic boundary con- 
ditions may introduce distortions in angular properties 
near the boundaries, since — as it will be shown in the 
following — some of the molecular correlators present un- 
expected peaks at wave vectors smaller than the center- 
of-mass peak position, corresponding to a distance of 
about three molecular diameters, and since the ampli- 
tude of these peaks grows at low temperature, we have 
made the additional effort to equilibrate a system made 
of 1728 molecules at T = 207K to make sure that no 
finite size effects show up in the molecular correlators for 
the system size considered. 

In order to calculate the Euler angles of each molecule 
the molecular reference frame has been chosen in such a 
way that the z-axes of the body-fixed frame has the same 
direction as the dipole of the molecule and the y-axes 
lies along the line connecting the two hydrogens. The 
x-axes is therefore determined to be orthogonal to the 
molecular plane. With this choice of the molecular axes 
the structure factors eq.(||) reduce, in the limit \q\ — > 
and for n, = n' = 0, to 



(31) 



(30) 



are the orientational correlation functions of the dipolar 
moments j2i(t) = //efi(i). In eq.(|3l|) Pi is the Legendre 
polynomial of order I and /j, is the dipole strength. These 
g-independent rotational correlation functions can be ex- 
perimentally measured for some values of I (see above). 
A complete analysis of the self part of C\ (t) for the sys- 
tem under investigation is treated in ref. Q). In what 
follows the discussion will be extended to the generalized 
correlators SJ^ ;/„/(<?, t) in the g-frame. 

In the first subsection we will discuss the static corre- 
lators and in the second subsection we will turn to the 
detailed analysis of the dynamic behavior of the corre- 
lation functions defined in eq.(|J) as calculated from the 
MD data. We will discuss the numerical results compar- 
ing them to the qualitative and quantitative predictions 
of MCT or MMCT. 



A. Static properties 

We have calculated all the static structure factors as 
defined in eq.(||) up to / = 2. The q- vector range has 
been chosen in such a way to include all the meaningful 
structure of each correlation function, i.e. the value of q 
ranges from 3.3nm _1 , which is the lower bound imposed 
by the finite size, up to around HOnm -1 , a value at 
which all the structure factors have practically reached 
their asymptotic values. The grid spacing has been fixed 
to Aq — l.llnm -1 which allows all the peaks in the 
correlators to be well resolved. 

In order to obtain a more significant statistics, the 
static structure factors have been calculated, for each 
configuration, for several directions of the q- vectors with 
respect to the fixed laboratory frame. The different re- 
sults are averaged after performing a suitable rotation 
which brings the g-vector along the z-axes of the lab 
frame (eq.(Q)). Furthermore the resulting correlation 
functions have been averaged over several configurations 
extracted from a time evolution which extends up to 
250ns. The global symmetries described in section II], 
which are all fulfilled within the numerical error by the 
numerical correlators, allow a further average of the static 
correlators. In this way the numerical noise is reduced to 
its lowest possible value. The resulting structure factors 
are shown in Figs, [j] and|| for the lowest simulated tem- 
perature T — 207 K. We found that within the numerical 
error, all correlators are real. 

Some of the static correlators shown in Figs, [j] and || 
are characterized by large amplitude peaks at q- vectors 
corresponding to distances of the order of three molecu- 
lar diameters. This is exactly the distance where finite 
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size effects could be expected for a system composed of 
6 3 molecules. To check if the results are real or simply 
artifacts of the simulation conditions, we report in Fig. 
H the static structure factors for the 6 3 and for the 12 3 
molecule system. We perform the check at the lowest 
studied temperature, where the angular correlations are 
enhanced. Although the statistic for the large system 
is poorer, as expected from the shorter simulation time, 
both position and amplitude of the peaks are not affected. 
Data in Fig confirms that the angular correlations do 
persist longer than the center of mass correlations and 
that the static molecular correlators, i.e. the quantities 
which are requested as input by the molecular mode cou- 
pling theory, are not affected by the size of the studied 
system. 

The different figures list the static correlation functions 
in terms of progressive angular complexity. Fig.0(a,b) 
includes all the diagonal correlators (i.e. I — I') with 
n = n' = 0. In Fig. 0(c) n and n' are still fixed to zero 
but the diagonality in I and I' is now relaxed. In Fig.0 
the correlators with n and/or n' different from zero are 
reported. Figs. [I] and show that several molecular cor- 
relators are comparable in intensity to the center of mass 
correlation function. Furthermore, some of them present 
completely new features, i.e. the q-dependence of the 
molecular correlations with I and I' different from zero 
is completely different from the center of mass structure 
factor. Some of the S]^ Vn ,{q) display peaks and min- 
ima at (j- vector values where the center of mass structure 
factor is rather structureless. The generalized spherical 
harmonics D l mn (flj) in eqs.(0) and (0) can be viewed as 
weights that "single out" molecules with certain orien- 
tations ttj. The differences in the g-dependence of the 
generalized orientational correlators show that looking at 
molecules with "selected" orientation reveals character- 
istic length scales of the system which are different from 
the center of mass ones. This can result in the shift of 
a peak as in the case of the main peak of S^q 10 (q) com- 
pared to SQ aQ (q) (see Fig. 0(a)). The "proper choice" 
of orientations can also reveal order on scales that are 
longer than the typical intermolccular distance and that 
is not visible in the center of mass correlators. For exam- 
ple, the most significant peak in Sl 20(9) an d Sl w (q) 
is located around q = 6.5nm _1 , a q- vector much smaller 
than the center of mass structure factor first peak. The 
q = 6.5nm~ 1 peak, which appears in correlators with 
/ or V equal to one and n = n' = 0, may reflects the 
strong dipolar interactions characteristic of water. At 
distances of the order of 27r/6.5 nm, about three molec- 
ular diameters, the SPC/E potential is equivalent to the 
potential generated by an electrostatic dipole with the 
same dipole moment of the SPC/E model. Thus, data in 
Fig. 0, suggest that in water angular correlations per- 
sist over distances much longer than the center of mass 
correlations. 
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FIG. 1. Static structure factors S;™ ;/„< (?)■ Diagonal cor- 
relation functions with n = n' = are shown in (a) and (b). 
Off-diagonal correlation functions with n = n' — are shown 
in (c). 

We also point out that the correlators in Fig. 0(c) have 
an intensity more than three times smaller than the ones 
reported in Fig. 0(a) (where for the latter one has to 
use as intensity S% l0 (q) — I). This suggests that the 
diagonal static correlation function could provide a good 
starting approximation for a MMCT description of the 
slowing down of dynamics in SPC/E water. We recall 
that this is at odd with the case of a liquid of Lennard 
Jones dumbbells |^6| where the off-diagonal terms have 
a large amplitude. Correlators with n different from zero 
(see Fig. 0) carry information about the planar shape of 
the molecule (n = n' = is equivalent to assuming that 
the Euler angle \ is always zero), i.e. on the absence of 
cylindrical symmetry. 
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FIG. 2. Static structure factors 5;™ ; / n /(g) with n and/or 
n 7^ 0. Subfigure (a) demonstrates the validity of eq.(p9|).The 
three upper lines refer to the case n = n' while the lower ones 
represent n = —n. 

Subfigures (a) of that figure demonstrates the validity 
of the property (|29|). Data in Fig. |^ shows that a few of 
these correlators describe a significant amount of angular 
correlations. Again, the largest amplitude is observed for 
the case in which I or V is one, i.e. for S\ 22(1) ■ 

Figure |2| also exhibits the n — > — n (or n! — > — n') sym- 
metry as given by eq.(p9[). Although this symmetry be- 
comes exact at least at high temperatures (see Appendix 
|a|), it seems to be valid within the numerical error even 
at T = 2Q7K, which is not high anymore. Although in 
Appendix [X] it is shown that for any temperature contri- 
butions to the distinct part exist for which this symmetry 
is still exact, it is not obvious to us why their weight is 
so large. 
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FIG. 3. Comparison of the static structure factors 
S^, v (?) at T = 207K between the 6 3 (solid lines) and the 
12 3 (symbols) molecules system. The comparison is presented 
for several S,™ Vn , (q), which have been arbitrarily shifted 
along the y-axis to improve the quality of the figure. From top 

to bottomiSio^o, 5^2, 22i 500,221 5 , 10 ,22i <5'l0,22j 5*22,221 5*22,221 

520,22 1 520,22 1 520,22 ■ The molecular correlators which have a 
peak at small q vectors are shown, to highlight the absence of 
finite size effects even in the worst case where the molecular 
correlation extends over about three molecular diameters. 

Figs.||, ||, H and show the temperature dependence of 
the static structure factors Vn , (q). In order to analyze 
the relative variation of the correlators on varying T we 
report three different temperatures, i.e. T = 207 K, 
T = 238 K and T = 400 K. Figs. | and | report the 
diagonal S"Q l0 (q). Fig. || reports the off diagonal terms 
with n = n' = 0, while Fig. |?] reports some correlation 
functions with ra^O. 
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FIG. 4. Temperature dependence of the diagonal static 
structure factors with I — 0, 1 and n = n' = 0. 



From these figures we see that the behavior of 
^ini'n'il) 011 changing T is strongly correlator depen- 
dent. The static structure factor Sq 0Q (q) (center of 
mass) shows a significant increase of the resolution of 
the peaks, which become sharper as T is lowered, espe- 
cially at small (/-vectors. The diagonal correlators with 

1 = I' = I are less temperature dependent. They remain 
almost unchanged at low temperatures changing T by 30 
degrees, from T — 238 K to T — 207 K, i.e. in the region 
where the molecular diffusivity decreases by more than 

2 order of magnitude |3jJ. Thus, while these correlators 
have a large amplitude, which implies that they might 
contribute substantially to the MMCT vertices, the dy- 
namical transition may not be controlled by them. 
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FIG. 5. Temperature dependence of the diagonal static 
structure factor with 1 = 2 and n = n' = 0. 

The off-diagonal terms shown in Fig. |^ as well as the 
static correlation functions with n and/or n' ^ (Fig. 0) 
also do not show significant temperature variation at low 
temperatures. This suggests that a small set of diago- 
nal correlators may play the relevant role in the slowing 
down of the molecular dynamics on supercooling. Of 
course, only a full MMCT calculation can confirm such 
hypothesis and can provide definitive answers on which 
modifications of the angular correlations drive the ideal 
glass transition. 
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B. Dynamic properties 

In this section we discuss the time evolution of the 
molecular structure factors with 1,1' < 2 at different 
g-vector values. We present the results of the calcula- 
tion in the perspective of testing how far the mode cou- 
pling framework can describe the behavior of S]™ Vn , (g, t). 
Since the complete MMCT equations have not yet been 
solved for the complete dynamic evolution of the correla- 
tors, we mostly keep our discussion to a qualitative level. 
We show that most of the universal predictions (see In- 
troduction) for the slow relaxation in supercooled liquids 
are excellently satisfied by the numerical data. 

One of the strongest predictions of MCT is the valid- 
ity of the so called time-temperature superposition prin- 
ciple in the a-relaxation region (see Introduction). The 
time-temperature superposition principle states that in a 
wide range of temperatures above the MCT critical point 
(which will be defined in the following) it is possible to 
scale the same correlator evaluated at different T on a 
single master curve through a rescaling of the time, i.e. 

<P(t) = 4>{t/r{T)) (32) 

In eq.(||^) cj>(t) indicates any dynamic structure fac- 
tor, <j)(t) the master function and r(T) is a temperature 
dependent time scale which is characteristic of the cho- 
sen correlator. The temperature dependence of the time 
scale is also predicted by MCT or MMCT. In leading or- 
der t(T) is a power law diverging at the critical transition 
temperature: 

t(T)~|T-T c |-t (33) 

Eq.(p3[) is one of the possible operative definitions of 
the mode-coupling critical temperature T c . T c is the tem- 
perature at which the characteristic time scales diverge, 
i.e the point at which the dynamic of the liquid is com- 
pletely frozen. Thus T c defines a kinetic transition from 
an ergodic to a nonergodic dynamic. This is a purely 
kinetic transition which does not have a thermodynamic 
counterpart. In real liquids no sharp transition is ob- 
served and close to T c the system switches to a different 
dynamic where hopping phenomena become dominant. 
Thus T c assumes the meaning of a crossover temperature. 
The transition temperature T c and the scaling exponent 
7 are not universal quantities since they strongly depend 
on the physical system under investigation and on the 
volume and pressure conditions, but they have a sort of 
"universality" in the sense that they are predicted to be 
correlator independent. Indeed, T c and 7 are predicted 
to have the same value for all correlation functions which 
couple to each other. In case of molecular liquids their 
values can be obtained from MMCT. In previous papers 
p3| , p4|j44| it has been shown that in the case of SPC/E, 
the time-temperature superposition principle is satisfied 
by self and collective center of mass correlators and also 



by the g-independent rotational correlation functions. It 
has been also shown that T c = 200 ±3 and the critical ex- 
ponent 7 = 2.7 are the same for all examined correlators 
within the numerical error. Here we generalize this con- 
clusion showing that the time-temperature superposition 
principle also holds for all the molecular collective corre- 
lators up to I = 2. We show here only two representative 
correlators. In Fig. || we have reported Si 0W (q, t) at 
the different temperatures normalized to its static value 
Sl ow (q). We have rescaled each curve choosing as t(T) 
the time at which the correlator has decayed to the value 
1 /e. The different curves perfectly overlap in the a-region 
confirming the validity of the time-temperature superpo- 
sition principle. The small graph shows as a reference 
the static structure factor; an arrow is pointing to the 
g-value for which the analysis is performed. In Fig. ^ 
we show the test of the time-temperature superposition 
principle for the same correlator but for a different value 
of q, while Fig. |l^ shows the same analysis for a different 
correlator, Sl io(q,t). These results, and similar analy- 
sis for the other correlators (not reported in this Article) 
lead to the conclusion that the time-temperature super- 
position principle is satisfied for all examined SjJJ ;/„/(<?)■ 

In order to verify the validity of the scaling law for the 
time scale r(T) and the "universality" of the exponent 
7 we represent in Fig. [ll] t -1 / 7 , for 7 = 2.7, as a func- 
tion of T. The two figures are at two different values of 
the q vector, i.e. q = l&nm^ 1 and q — 28nm~ 1 and in 
each figure t _1 / 7 (T) is shown for all examined correla- 
tion functions (i.e. I < 2, — I < m < I and —l<n< I). 
Data in Fig. [ll] suggest that, as previously observed for 
the center of mass and g-independent rotational correla- 
tors, the power law ( p3| ) is well satisfied with the same 
"universal" values of the critical exponent 7 = 2.7 and of 
the critical temperature T c = 200 ± 3. We also note that 
the time scales for fixed T vary by about one decade as 
can be observed by the rather different slopes. This anal- 
ysis strongly confirms the MCT and MMCT prediction of 
the existence of a unique critical temperature T c at which 
both translational and rotational degrees of freedom cross 
from an ergodic to a nonergodic dynamics following a 
power law behavior ruled by an "universal" exponent 7. 
Thus, MMCT seems to be a good framework in order to 
describe qualitatively and semi-quantitatively the tem- 
perature dependence of the structural correlations in a 
supercooled molecular liquid in a wide range of temper- 
atures above a critical T c , as defined by eq. (|3^) . 

Data in Fig. |j, [)] and [n] show oscillations at short 
time connected with the librational and vibrational dy- 
namics of the water molecules. These motions (with a 
characteristic timescale of the order of 0.1 ps) modulate 
the approach to the plateau and may interfere with the 
universal dynamics characteristic of the /9-region p5j . 
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FIG. 11. a-relaxation time to the power — 1/7 for all correlation functions with / < 2, — I < m < I and — / < n < I at two 
different g-vector values. Lines are drawn to guide the eyes. 
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If the temperature is lowered very close to T c , MCT 
provides quantitative predictions for the time evolution of 
the correlation function. We stress again that these pre- 
dictions are the same for MMCT. In the region of the first 
scaling law MMCT predicts for 5^ ; , n ,(g, t) the factor- 
ization of the time-dependence from the dependence on 
q, I, n, I', n! and to, i.e. near to T c its time-dependence is 
given by the so-called /3-correlator G(t) which describes 
the dynamics close to the plateau on a time scale t a (T) 
(/3-regime). The equation for G(t) does not depend on 
y, I, n, I', n' and to, it only involves the exponent param- 
eter A. This equation which is the same for MCT and 
MMCT can be solved exactly in the asymptotic limits 
t <C t a and t 3> t a yielding for both cases a power law 
dependence with exponent a (critical law) and exponent 
b (von Schweidler law), respectively. The first one de- 
scribes the relaxation onto and the second one from the 
plateau. 

The two exponents are both functions of A according 

to M 



A = 



r(i -a) 2 _ r{i + b) 2 



(34) 



r(l-2a) r(l + 26) 

where T is the Euler gamma function, and they are con- 
nected to the relaxation time exponent 7 by the relation 

1 1 

Ta + 2b 



MMCT provides an explicit expression for A which con- 
tains the static molecular correlators at T c p8| . 

The von Schweidler law which also describes the early 
a-relaxation regime is given by 

(36) 



-O(P) 



7 



(35) 



where F™'f, n ,(q), ff^™ , (q) and H$™ n , (q) are, respec- 
tively, the critical nonergodicity parameter, the critical 
amplitude and the amplitude of the next to leading order 
correction for the unnormalized correlators. An explicit 
expression for h^ 2 '(q) = H^ 2 \q) / S(q) for simple liquids 
has been recently derived |47|| . The time window in which 
the asymptotic power law t b holds is strongly correlator 
and q dependent. Hence to perform a careful MCT anal- 
ysis it is always necessary to take into account also the 
second order corrections. This has previously been found 
for the center of mass correlator for water p3 and for the 
molecular correlators for diatomic molecules po| . In or- 
der to test the validity of the von Schweidler law and the 
relevance of the second order corrections we have fitted 
the time evolution of our correlators according to eq. ( |3^ ) 
for the lowest simulated temperature T — 207K which is 
few degrees above the critical one. 
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In Figs. |12|- |15| we show the fits with eq.(|3q) performed 
for some representative normalized diagonal correlators. 
The vertical lines indicate the time window selected for 
the fit which has been chosen consistently with the pre- 
vious fits for translational and (/-independent correlators, 
i.e t covers a range of two orders of magnitude, from 7ps 
up to 800ps. Fig. ^represents the fit with eq.(|6j) using 
the value b = 0.5 as previously found in the MCT study of 
the time dependence of self and collective center of mass 
correlators. We show the correlator Sj 10 (g,t) for three 



different values of the q vector, marked by arrows in the 
small graph representing the static Sl 010 (q). The solid 
line is the numerical curve, while the long dashed line 
is the result of the fit. The quality of the fit is remark- 
ably good. The g-dependence of the validity of the fit is 
evident noting that for q = 8rim _1 and q = 15.5nm _1 
the fitted von Schwcidler law holds far above the fitting 
range, while for q = 36.5nm _1 the two curves separate 
out of the fit window. 
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To clarify the importance of the second order correc- 
tions we have also reported in the figure the curve ob- 
tained keeping only the term t (short dashed line). We 
see that while in the case q = 15.5nm _1 the first power 
law fit alone extends to almost three order of magnitude, 
for the other values of q the second order corrections are 
necessary. The obtained fitting parameters F, and 
are shown in Fig. [l6] as functions of q together 
with the static structure factors. Analytic expressions 
for the calculation of these quantities are provided by 
MMCT (for the amplitude iJ (2) no explicit expression 
for the molecular system has been derived yet), so that 
they can be, in principle, calculated using the molecular 
static structure factors. For diatomic molecules this has 
been done for F and ijM §§]. 

The fit to other correlators show similar behaviors. We 
want to stress that, according to MCT and MMCT pre- 
dictions, in all the fits to eq.(|3^) the same exponent b 
has turned out to be satisfactory. The independence of 
the power law exponent b on the chosen correlator and 



on the q value is a strong argument in favor of MMCT 
as framework to describe the slow relaxation in super- 
cooled molecular liquids. We also stress that the two 
independently calculated exponents b and 7 satisfy the 
theoretical prediction of eq.(|35|). 

In view of a future comparison with the full time de- 
pendence of the MMCT correlations, we report in Fig. [l?] 
the parameters of the fit to the numerical correlation 
function according to a stretched exponential form 



exp 



(37) 



We note that in all examined cases the large g-vector 
limit of the stretching parameter /3 is about 0.5, i.e. equal 
to the value of b. Such equivalence is predicted by MCT 
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V. SUMMARY AND CONCLUSIONS 
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In this paper we have analyzed the slow relax- 
ation properties of a supercooled liquid of planar rigid 
molecules described by means of a MD simulation. The 
parameters of the simulation have been chosen in such a 
way as to mimic the thermodynamic and dynamic prop- 
erties of liquid water, a molecular liquid characterized by 
a network structure at low temperatures. Starting from 
the set of MD data we have calculated all the static and 
dynamic molecular structure factors S{£ Vn , (q, t) up to 
I = 2 in the q-frame. These correlators present a va- 
riety of features which open a deeper insight into the 
structural properties of the liquid. For example the sym- 
metries of the rotational correlators reflect the geomet- 
rical properties of the water molecules, while the peak 
distribution describes the intermolecular interactions. It 
becomes clear from an overview on the complete set of 
correlators that a theory able to describe the relaxational 
properties of a supercooled molecular liquid can not ne- 
glect the molecular correlators with 1,1' ^ and their 
coupling to the center of mass one. For the static dis- 
tinct part at least at high temperatures we have proven 
that it is independent of the sign of n and n' for molecules 
with C rv symmetry. The MD-results indicate that the 
independence even holds at lower temperatures (down to 
207K) within the numerical error. 

In this Article we argue that MCT in its molecular 
formulation is a good candidate to be the correct frame- 
work for the description of the slow dynamics in super- 
cooled molecular liquids. In order to support this view 
we have compared the properties of our MD liquid to 
the universal predictions of MCT (or MMCT) which 
are independent of the molecular nature of the liquid. 
Even if the complete set of MMCT equations has not yet 
been solved, thus preventing from a rigorous comparison 
between theory and simulation, we have tested several 
asymptotic MCT predictions. We have shown that the 
temperature and time dependence of our MD correlators 
perfectly obeys the time-temperature superposition prin- 
ciple predicted by MCT and MMCT, and we have found 
the MMCT "universality" , i.e. the critical temperature 
T c the exponents b and 7 (and therefore also a) do not 
depend on q, I, n, V, n' and m. 

Close to T c , where the asymptotic /3-correlator behav- 
ior is reinforced, we have compared the numerical cor- 
relator with the predicted von Schweidler law plus sec- 
ond order corrections in the early a-region. We have 
shown good quality fits consistent with what was previ- 
ously found for the self and collective dynamics of the 
center of mass alone and of the ^-independent angular 
correlators. Again, the predicted relation ( |35| ) between 
the von Schweidler exponent and the relaxation time ex- 
ponent is fulfilled, as well as the large g-vector limit of 
all stretching exponent /3 ; ™ Vn ,{q)- 
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APPENDIX A: MOLECULES WITH 
Crv -SYMMETRY 

In this appendix we will investigate the validity of the 
n — » — ri symmetry for the static distinct part of the cor- 
relators (c.f. eq.(p6|)) for molecules with C r „ -symmetry. 

We consider a liquid of N identical, rigid molecules. 
Each molecule is made up of N a atoms of type a, 
a = 1,2,.... The position of atom v a of type a in 



the jth molecule is denoted by 2y„ a , j 
v a = 1,2, ...,N a . Then, in a site-site representation the 
total potential energy V is given by: 



1,2, 



,N, 



(Al) 



where we restrict ourselves to two-body interactions with 
pair potential v aa i(x) between atoms of type a and a'. 
Introducing the microscopic molecular density for atoms 
of type a in molecule j: 



N,-, 



(a) 
Pj 



and its Fourier transform: 



(a) 



(A2) 



(A3) 



!/„ = ! 



we can rewrite (Al) as follows: 
1 



(x)v aa > (x - x')p)? ' (x') 



2 EE 



d 6 x I d 3 x' 



(<*')< 



(A4) 



or by use of dA3| ): 



2 tri , v 

j=F3 q 



(A5) 



with v aa '(q) = J d 3 x v aa '(x)e lq ' x 1 the Fourier transform 
of the pair potential which depends on q — \q\ and the 
volume V of the system only. 
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Instead of the site-site coordinates we introduce center 
of mass and relative coordinates: 



and 



Ja) _ J.a) _ -j 



(A6) 



(A7) 



respectively. m a i s the mas s of a tom s of type a. Next 
we substitute (A6) and (A7) into ( [A3] ) which leads to 



where 



(A8) 



(A9) 



Now, let } be the relative vectors in the body hxed 
frame and flj = (<j>j, 9j, Xj) the Euler angles between the 
laboratory and the body fixed frame. Then it is: 

^l=R- l {^)R^ (Aio) 



with R(fl) G 50(3). Substitution of ( |A1O0 into flAg) and 
making use of the Rayleigh expansion |pq| we arrive at 



N a 



£ yu%2) 



(All) 



and with the transformation of the spherical harmonics 
under R G SO (3) ||: 



III) 



N a 



£ *w4? 



(A12) 



ii(r) are the spherical Bessel functions, Dj n , m (f2) are 
Wigner's rotation matrices, R a — l-Rw" | for all u a , 
eh*) = Ri^/Ra and e q = q/q. The substitution of 



pf\q) @ with pf\q,Clj) from JIT^) into (jAf) yields 
the center of mass-angular representation of V . 

V{{af"l ' }) = V(£ 1> ... 1 2 Jv; n 1> ...,fi JV ) = 



9£t/£ £ Vlmn,l> m >n'{q) ' 



. e -^(^-»^ n (^.)^(%0 (A13) 
with the transformed pair potential 



Vl mn ,l'm'n'(q) = (47r) 2 ^ 'l^ (e g ) Y{! m , (e q ) }^ cdot 

a a.' 

, {q)ji{qRa)ii i {qRa') ■ 

i ( A i4) 



£^(et ) )] ( £ 1W(# 



The representation of y by ( Al3| ) is quite obvious, 
since the interaction potential Vq(x, 0, fi') between iwo 
molecules with center of mass separation x and orien- 
tation fl and 17' can be expanded with respect to the 
complete set of functions e^ s D l mn (fl)D l m * n , (ft'). 

All the local molecular symmetry is contained in the 
coefficients Vi mn> i' m ' n '(q) through the entity: 



Hi 



(a) 



£ ^n(4: 3 ). 



(A15) 



which will be discussed for C rv symmetry. In the follow- 
ing we choose the body fixed z-axis along the molecular 
r— fold symmetry axes. First of all we note that the r-fold 
rotational symmetry implies: 



•'// 



(a) 



^ , n G {0,±r,±2r, ...} 
= , otherwise 



(A16) 



Here we have assumed that all types of atoms a = 1, 2, ... 
have exactly an r-fold rotational symmetry. This is not 
necessary for a molecule with r-fold axes. There may 
be some types of atoms which have a 2r-, 3r— , etc. fold 
symmetry. For such a it is j/j" for n = p^r where 



2,3,.... Nevertheless, ( A16| ) remains true since 



p = min Q p^ — 1, by assumption. 

Now let us turn to the reflection symmetry. With 
(4>)% , (where #( Q ) does noi depend on i/ Q , due to 
the C rv symmetry and the appropriate choice of the body 
fixed z-axes) we denote the azimuthal and polar angles 
of ei" 3 . For r even and fixed a the angles can be 
chosen such that the reflection symmetry with respect to 
the (a; — z)-plane implies: 



7T 



(a) 



A„ 2tt 



(A17) 



for v a = 1, 2, JV a . The integer A„ a may depend on 
i/„ . using the explicit (0, 0) dependence of Y/ n (e) = 



H„(0,$ = c hl P ln (0)e m * §§, we get from Q with 



(A17) 



ClnPln(0 {a} )(-l) n £ e" 

(-i)"q»P/„(0 (q) ) 2 e ~ 



■m0 



■in<fii, 
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N„ 



= E YiM:\o {a} ) 



-yt ] (Ais) 

-n. In the case that r is odd a similar choice 



where n ■ 

o{(f>^} can be done, such that (A18) holds. Due to (AIS) 



we find that the coefficients vimn,i'm'n' (q) do not depend 
on the sign of n and n' . Therefore (A13) can also be 
rewritten as follows: 

V(xi, ...,x N ;Qi, ...,Qn) = 



1 



1 



q kk' 

+ £>„(%)] [D K ,(Sl f ) + D^,(n r )}* (A19) 



3+3 9 KK 



where k = (l,m,n), k = (Z,m, n) and D K = -D^ m has 
been used as a shorthand notation. ( Al9| ) is the basic 
result which allows to study the n — > — n symm etry for 
the static distinct part of the correlators. From (A6) we 
get for t = 0: 

■SgUfl = i''"' t(2i + l)(2r + 1)]' ^ E ^ ■ 

AT 

. e -/3V(iEi,...,:Ejv;ni,...,f2jv) (A20) 

where Z c is the configurational partition function. Now 
we expand exp(— f3V) into a power series. The zeroth 
order term of Sjf K ,(q) is proportional to 8 n ^5 n >Q which in 
a trivial sense is invariant under n or n' into —n or —n', 
for arbitrary q,l,m,l' and m'. For the first order term 



we find with flA19p : 



order 



(A21) 



where the factor of proportionality is independent of 
to, to', n and n'. Since u KjK ' does not depend on the sign 
of n and n' for all q, Z, to, Z' and to', the same is true for 
S^fl,(q) in first order. In second and higher order prod- 
ucts of the rotation matrices occur. In the /u-th order 
contribution of type 



1 

r 7 - 



E E 

91 — 81/ "l—"'' 



.(&) 



• IK ;<),)•..../),. ;<>,: ■£>;,(%0.....^(^.,) (A22) 

occur with v < \i. For f = 1 it coincides with the first 
order contribution. Using the product rule for the -D's, 
we get for v — 2 a contribution to S^ K , (q) which is pro- 
portional to: 



A kI>($) = yE E C(Kl«2«) C («l«2«') ■ 



9192 «i> 



(?l)U«2«i (92) 



(A23) 



where Y^ qi q 2 denotes summation over q\ and qi such that 
9i + 92 = q and 

C(kiK2k) = C{l\lil\ m\m2rri)C{l\lil'-, rixn^n) (A24) 

is a product of Clebsch Gordan coefficients. Then it fol- 



lows from (A.23) with the properties of C(hl2l; nin 2 n) 



A fAo) = ^E' E (-i) ,1+,!+, c(^2 K )cK4«') 



^k;®)^^ (92) 



(A25) 



where we have used the independence of v KlK / and of 
Vk 2 k' 2 0n the sign of n\ and First of all we observe 
that all terms in the sum of ( A25) with Z1+Z2+Z even coin- 
cide with the corresponding terms in (A22), independent 
of q, toi and to-2- Hence for those terms the n — > — n 
symmetry holds for arbitrary q,h,mi,l2 and TO2, pro- 
vided Zi + Z2 + Z is even. On the other hand choosing 
to = and q = qo = (0, 0, (7) (g-frame) and changing 
toi — » —mi, TO2 — > — TO2 we find 

4«'® = t?E' E c( Kl «2K)c(44«') ■ 



9192 



(A26) 



where k — (l,m,n). The (gxTOi, 92^2)— dependence of 
(51)^2^ (92) is given by (cf. ( |A14| ) ): 



1 l 1 m 1 V-91 y - 1 (2«i 2 V c 92 / ^ L 



(A27) 



where = to = toi + TO2 has been used. The azimuthal 
angles 4>\ and (f>2 of e qi and e 92 fulfill cf>2 — 4>i = 7r due to 
9oe qi + 9oe q2 = 0, i.e. the r.h.s. of ( A27| ) equals (— l) mi 
which does not depend on the sign of mi. Therefore 
we find that Yi lSll (e 9l )^ 2 m 2 (e 92 ) = ij imi (e 9l )Y hm2 (e q2 ) 
and accordingly that 



w«i<(«!)% 2 «i(92) = «ki<(9i)^«2« 2 (92)- 
Making use of (|A28|) , we get from (|A26|) that 

4« ) '(?o)=4 2 2'(?o) 



(A28) 
(A29) 

for all k' and k with to = 0. We see that the n — > — n (or 
n' — > — n') symmetry of all the second order contributions 
only holds in the g-frame and for to = to' = 0. 
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(2) 

Besides A\.^, (q) , in third order there also exist contri- 
butions of type 



4^ = ^ £' 



E 



(&) (A30) 

where (...) contains products of C's, and the sum over 
<fi , 52 an d g 3 is restricted to <fi + 92 + <73 = 5- For q = Qq 
and m = the same procedure as for v — 2 can be used 
in order to prove that the n — > —n symmetry will be true 
if the three-fold product of the v's in ( A29j ) is invariant 
under m, — » -m„ i = 1,2,3. This requires that (in 
analogy to v = 2): 



-i(m 1 cf> 1 +m2'p2+m a <p a ) _ i(m 1 (j> 1 +m 2 cl>2+m 3 <p 3 ) 



(A31) 



Although = m = mi + 1x12 + m^, the azimuthal angles 
of e qi (in contrast to v = 2) are not restricted because 
Si=i <7o • % = does not fix ^3 if <fii and 02 are given. 
Therefore the generic case for q\ obeying 53i=i <7o ' <fi = 
will not satisfy (A30). But similar to v = 2, there 



is a subset of terms in the sum of (A29) which will be 
invariant under n — > —n (or n' — > — n') for a/Z q,l,m,l' 
and m'. 



The results show that the invariance of S, 



(d) 

Imn.l'm 



for arbitrary (q,l,m,l' ,m') and for (<f = qb,l,l',m — 
m! = 0) only holds up to, respectively, the first and sec- 
ond order in 1/ksT. 
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